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1. - Introduction: Localized solutions to the wave equations 

Since 1915 Bateman[l] showed that Maxwell equations admit (besides of the ordinary 
planewave solutions, endowed in vacuum with speed c) of wavelet-type solutions, en- 
dowed in vacuum with group-velocities < v < c. But Bateman's work went practically 
unnoticed. Only few authors, as Barut et al.[2] followed such a research line; inciden- 
tally, Barut et al. constructed even a wavelet-type solution travelling with Superluminal 
group-velocity[3] v > c. 

In recent times, however, many authors discussed the fact that all (homogeneous) 
wave equations admit solutions with < v < oo: see, e.g., Donnelly & Ziolkowski[4], 
Esposito[4], Vaz & Rodrigues[4]. Most of those authors confined themselves to investi- 
gate (sub- or Super-luminal) localized non-dispersive solutions in vacuum: namely, those 
solutions that were called "undistorted progressive waves" by Courant & Hilbert. Among 
localized solutions, the most interesting appeared to be the so-called "X-shaped" waves, 
which — predicted even by Special Relativity in its extended version[5] — had been mathe- 
matically constructed by Lu & Greenleaf[6] for acoustic waves, and by Ziolkowski et al.[7], 
and later Recami[8], for electromagnetism. 

Let us recall that such "X-shaped" locahzed solutions are Superluminal (i.e., travel 
with V > cin the vacuum) in the electromagnetic case; and are "super-sonic" (i.e., travel 
with a speed larger than the sound-speed in the medium) in the acoustic case. The first 
authors to produce X-shaped waves experimentally were Lu & Greenleaf[9] for acoustics, 
and Saari et al.[10] for optics. 

Notwithstanding all that work, still it is not yet well understood what solutions (let 
us now confine ourselves to Maxwell equations and to electromagnetic waves) have to 
enter into the play in many experiments. 

2. — About evanescent waves 

Most of the experimental results, actually, did not refer to the abovementioned local- 
ized, sub- or Super-luminal, solutions, which in vacuum are expected to propagate rigidly 
(or almost rigidly, when suitably truncated). The experiments most after fashion are, on 
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the contrary, those measuring the group-velocity of evanescent waves[cf., e.g., refs.11,12]. 
In fact, both Quantum Mechanics [13] and Special Relativity [5] had predicted tunnelling 
wavepackets (tunnelling photons too) and/or evanescent waves to be Super luminal. 

For instance, experiments[12] with evanescent waves travelling down an under- 
sized waveguide revealed that evanescent modes are endowed with Superluminal group- 
velocities[14]. 

A problem arises in connection with the experiment [15] with two "barriers" 1 and 
2 (i.e., segments of undersized waveguide). In fact, it has been found that for suitable 
frequency bands the wave coming out from barrier 1 goes on with practically infinite 
speed, crossing the intermediate normal-sized waveguide 3 in zero time. Even if this 
can be theoretically understood by looking at the relevant transfer function (see the 
computer simulations, based on Maxwell equations only, in refs.[16,17]), it is natural to 
wonder what are the solutions of Maxwell equations that can travel with Superluminal 
speed in a normal waveguide (where one normally meets ordinary propagating — and not 
evanescent — modes) . . . 

Namely, the dispersion relation in undersized guides is u;"^ — k"^ — — Q^, so that the 
standard formula v ~ dcu/dk yields a v > c group- velocity [17, 18]. However, in normal 
guides the dispersion relation becomes u;'^ — k'^ — +0,'^, so that the same formula yields 
values V < c only. 

We are going to show that actually localized solutions to Maxwell equations propa- 
gating with ^; > c do exist even in normal waveguides; but their group-velocity v cannot 
be given*^ by the approximate formula v ~ da;/dA;. One of the main motivations of 
the present note is just contributing to the clarification of this question. 

3. — About some localized solutions to Maxwell equations. 

Let us start by considering localized solutions to Maxwell equations in vacuum. A 

theorem by Lu et al.[19] showed how to start from a solution holding in the plane {x,y) 

for constructing a thrcedimensional solution rigidly moving along the z-axis with Supcr- 

Let us recall that the group- velocity is well defined only when the pulse has a clear bump in space; 
but it can be calculated by the approximate, elementary relation v ~ dw/dfc only when some extra 
conditions are satisfied (namely, when a; as a function of k is also clearly bumped). 
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luminal velocity v. Namely, let us assume that ip{p;t), with p = {x,y), is a solution of 
the 2-dimensional homogeneous wave equation: 

{9l + d'y-^df) i^{p-t) = 0. (1) 

By applying the transformation p ^ p sinO; t ^ t — {cos 9 /c)z, the angle 9 being 
fixed, with < ^ < 7r/2, one gets[19] that ■0(p sin^; t — (cos^/c) z) is a solution to the 
threedimensional homogeneous wave-equation 

(V^ - V'(p sin^; t-^z^^O. (2) 

The mentioned theorem holds for the vacuum case, and in general is not valid when 
introducing boundary conditions. However we discovered that, in the case of a bidimen- 
sional solution ip valid on a circular domain of the {x, y) plane, such that ^0 = for |p| =0, 
the transformation above leads us to a (three-dimensional) localized solution rigidly trav- 
elling with Superluminal speed v — c/ cos9 inside a cylindrical waveguide] even if the 
waveguide radius r will be no longer a, but r — a/ sm.9 > a. We can therefore obtain 
an undistorted Superluminal solution propagating down cylindrical (metallic) waveguides 
for each (2-dimensional) solution valid on a circular domain. Let us recall that, as well- 
known, any solution to the scalar wave equation corresponds to solutions of the (vectorial) 
Maxwell equations (cf., e.g., ref.[8] and refs. therein). 

For simplicity, let us put the origin O at the center of the circular domain C, and 
choose a 2-dimensional solution that be axially symmetric ^(p;t), with p = \p\, and 
with the initial conditions ■0(p;t = O) = 0(p), and dt/^/dt = ^{p) at t — 0. 

Notice that, because of the transformations 



p =^ p sin 9 (3a) 

COS^ /„! N 

t^t z , (3b) 

c 

the more the initial ip{p; t) is localized at t = 0, the more the (threedimensional) wave 
ip{p sin 9; t — {cos9/c)z will be localized around z = vt. It should be also emphasized 
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that, because of transformation (3b), the velocity c goes into the velocity v — c/cos9 > c. 
Let us start with the formal choice 



m = ^ ; ^ . (4) 

p 



In cylindrical coordinates the wave equation (1) becomes 

(^^d,pd, - la,^) i;{p;t) = 0, (!') 

which exhibits the assumed axial symmetry. Looking for factorized solutions of the type 
ip{p;t) = R{p)-T{t), one gets the equations d^T ^ -lo^T and {p'^dp+d'^+u;^ /c^)R ^ 0, 
where the "separation constant" a; is a real parameter, which yield the solutions 

T — A cos Lot + B sin Lot 

(5) 

R = C Jo(-p) , 
c 

where quantities A,B,C are real constants, and Jq is the ordinary zero-order Bessel 
function (we disregarded the analogous solution Yo{Lop/c) since it diverges for p — 0). 
Finally, by imposing the boundary condition = at p = a, one arrives at the base 
solutions 

k uj 
ilj{p;t) = Jo{—p) {Ancosoont + BnSinujnt) ; k=-a, (6) 

Gj c 

the roots of the Bessel function being 

_ ^nd 
Kjl — • 



The general solution for our bidimensional problem (with our boundary conditions) 
will therefore be the Fourier-type series 

k 

*2D(p;i) = E~=i^o(— p) (Ancosojj + B^ smooj) . (7) 
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The initial conditions (4) imply that X] A„ Jo(A;„p/a) = S{p)/p, and YjBnJo{knp/o) — 
0, so that all Bn must vanish, while An — 2[a?Jl{kn)]~^] and eventually one gets: 

where ujn — knc/ a. 

Let us explicitly notice that we can pass from such a formal solution to more physical 
ones, just by considering a finite number N of terms. In fact, each partial expansion will 
satisfy (besides the boundary condition) the second initial condition dt'4> = for i = 0, 
while the first initial condition gets the form = /(p), where f{p) will be a (well) 
localized function, but no longer a delta-type function. Actually, the "locahzation" of 
0(p) increases with increasing N. We shall come back to this point below. 

4. - Localized waves propagating Super luminally down (normal-sized) waveg- 
uides. 



We have now to apply transformations (3) to solution (8), in order to pass to 
threedimensional waves propagating along a cylindrical (metallic) waveguide with radius 
r — a/ sm.9. We obtain that Maxwell equations admit in such a case the solutions 

*3D(p,^;t) = E^=i ( ^^|r|-^ J Jo(— psin^) cos 

where ujn — knc/a, which are sums over different propagating modes. 

Such solutions propagate, down the waveguide, rigidly with Superluminal velocity*^ 

V = c/cos6. Therefore, (non-cvancscent) solutions to Maxwell equations exist, that are 

waves propagating undistortcd along normal waveguides with Superluminal speed (even 

if in normal-sized waveguides the dispersion relation for each mode, i.e. for each term of 

the Fourier-Bessel expansion, is the ordinary "subluminal" one, uJ^ j(? — k\ = +il^). 

It is interesting that our Superluminal solutions travel rigidly down the waveguide: 

Let us stress that each eq.(9) represents a multimodal (but localized) propagation, as if the geo- 
metric dispersion compensated for the multimodal dispersion. 



kn COS 9 , c , 

a cos 6 



(9) 
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this is at variance with what happens for truncated (Superluminal) solutions [7- 10], which 
travel almost rigidly only along their finite "field depth" and then abruptly decay. 

Finally, let us consider a finite number of terms in eq.(8), a,t t — 0. We made a few 
numerical evaluations: let us consider the results for N — 22 (however, similar results can 
be already obtained, e.g., for N — 10). The first initial condition of eq.(4), then, is no 
longer a delta function, but results to be the (bumped) bidimensional wave represented 
in Fig.l. 

The threedimensional wave, eq.(9), corresponding to it, i.e., with the same finite 
number N — 22 oi terms, is depicted in Fig.2. It is still an exact solution of the wave 
equation, for a metallic (normal-sized) waveguide with radius r — a/ sin^, propagating 
rigidly with Superluminal group- velocity v — c/ cos9; moreover, it is now a physical 
solution. In Fig.2 one can see its central portion, while in Fig. 3 it is shown the space 
profile along z, ior t — const., of such a propagating wave. 
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Figure Captions 



Fig.l — Shape of the bidimensional solution of the wave equation vahd on the circular 
domain p < a; a — 0.1 mm of the {x,y) plane, for t — 0, corresponding to the sum of 
N — 22 terms in the expansion (8). It is no longer a delta function, but it is still very 
well peaked. By choosing it as the initial condition, instead of the first one of eqs.(4), one 
gets the threedimensional wave depicted in Figs. 2 and 3. The normalization condition is 
such that |*2d(p = 0; i = 0)p = 1. 

Fig.2 - The (very well locahzed) threedimensional wave corresponding to the initial, 
bidimensional choice in Fig.l. It propagates rigidly (along the normal-sized circular 
waveguide with radius r — a/ sinO) with Superluminal speed v — c/ cos 9. Quantity 77 is 
defined as 77 = {z — ^^t). The normalization condition is such that |^'3d(p = 0; 77 = 
0)P = 1. 

Fig. 3 - The shape along z, ai t — 0, oi the threedimensional wave whose main peak is 
shown in Fig.2. 
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